Diffusion of gold nanoclusters on graphite 
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We present a detailed molecular-dynamics study of the diffusion and coalescence of large (249-atom) 
gold clusters on graphite surfaces. The diffusivity of monoclusters is found to be comparable to that 
for single adatoms. Likewise, and even more important, cluster dimers are also found to diffuse at 
a rate which is comparable to that for adatoms and monoclusters. As a consequence, large islands 
formed by cluster aggregation are also expected to be mobile. Using kinetic Monte Carlo simulations, 
and assuming a proper scaling law for the dependence on size of the diffusivity of large clusters, we 
find that islands consisting of as many as 100 monoclusters should exhibit significant mobility. This 
result has profound implications for the morphology of cluster-assembled materials. 

PACS numbers: 36.40. Sx, 61.46. +w, 68.35.Fx, 07.05. Tp 



I. INTRODUCTION 

Nanometer-size clusters — or simply nanoclusters — 
are intrinsically different from bulk materials. Elfl Yet, un- 
derstanding of several of their most fundamental physical 
properties is just beginning to emerge (see for instance 
Refs. ^-^), thanks largely to rapid progresses in the tech- 
nology of fabrication and analysis, but also considerable 
advances in computational tools and pjethodology. 

It has recently been demonstratedEElilj that deposit- 
ing clusters (rather than single atoms) on surfaces allow 
the fabrication of interesting nanostructured materials 
whose properties can be tailored to specific technolog- 
ical applications, e.gp-,micro-electronic, optoelectronic, 
and magnetic devices.Ej If single-atom deposition is used, 
the nanostructures have to be grown directly on the sub- 
strate through diffusion and agregation, which depends 
in a detailed (and in general very complicated) manner 
on the interactions between surface atoms and adatoms. 
By contrast, for cluster deposition, the clusters are pre- 
pared before they hit the surface, giving considerably 
more flexibilitytS in assembling or organizing clusters 
for particular applications. It has been shown, for in- 
stance, that by changing the mean size of the incident 
carbon clusters, it was possible to modify the structure 
of tlMi resulting carbon film from graphitic to diamond- 
like.t3 This however requires that sufficient control over 
the clustePndeposition and subsequent growth process be 
achieved.tllll3 

Diffusion evidently plays a central role in the fabrica- 
tion of thin films and self-organized structures by clus- 



ter deposition. It has been demonstrated experimentally 
that gold or antimony clusters diffuse on graphite sur- 
faces at a surprisingly high rate of about 10~^ cm^/s at 
room temperature|l3 quite comparable to the rates that 
can be achieved by single atoms in similar conditions. 
This was confirmed theoretically Jay Deltour et al. using 
molecular-dynamics simulations :L| clusters consisting of 
particles which are incommensurate with the substrate 
exhibit very rapid diffusion. The cluster diffuses "as a 
whole" , and its path is akin to a Brownian motion in- 
duced by the internal vibrations of the clusters and/or 
the vibrations of the substrate. This is in striking con- 
trast with other cluster diffusion mechanisms, whereby 
the motion results from a combination of single-atom 
processes, such as evaporation-condensation, edge diflFu- 
sion, etc. The latter mechanisms are more appropriate 
to clusters which are in epitaxy with the surface, and are 
likely not significant in cases where the mismatch is large 
and/or the substrate-cluster interactions are weak, such 
as in Refs. ^ and [l^ (see ReL ^ for a review). 

In the present paper, we re-examine the problem 
of cluster diffusion in the cluster-substrate-mismatched 
case, now using a much more .-accurate model: indeed, 
in the work of Deltour et al.p cluster-cluster, cluster- 
substrate and substrate-substrate interactions were all 
assumed to be of the Lennard-Jones form, which can- 
not be expected to correctly describe "real materials". 
Here, we consider a simple, but realistic model for the 
diffusion of gold clusters on a graphite surface (HOPG). 
We are concerned with gold becausadtias^ been the ob- 
ject of several experimental studies,E30il3 but also be- 
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cause realistic semi-empirical, many-body potentials are 
available for this material. The energetics of gold atoms 
is desciibed in terms of the embedded-atom-method 
(E AM), E3 while c asbon atoms are assumed to interact via 
Tersoff potentials the (weak) interactions between gold 
and carbon atoms are modeled with a simple Lennard- 
Jones potential. A comparable model was used recently 
by Luedtke and Landman to study the |anomalous dif- 
fusion of a gold nanocluster on graphitep diffusion was 
found to proceed via a stick-slip mechanism, resulting in 
an apparent Levy-flight type of motion. In the present 
work, we examine closely the variations with tempera- 
ture of the rate of diffusion, as well as the microscopies 
of cluster dimers (diclusters). 

We find the diffusivity of monoclusters to be entirely 
comparable to that for single adatoms. Likewise, and 
most important, diclusters are also found to diffuse at a 
rate which is comparable to that for adatoms and mon- 
oclusters. It is therefore expected that large islands, 
formed by the aggregation of many clusters, should also 
be mobile. Based on this observation, we carried out 
kinetic Monte Carlo simulations of island diffusion and 
coalescence assuming a proper scaling law for the depen- 
dence on size of the diffusivity of large clusters. We find 
that islands consisting of as many as 100 monoclusters 
exhibit significant mobility; this is consistent with the 
observation on graphite of large (200 monoclusters) gold 
islands. The morphology of cluster-assembled materials 
is profoundly affected by the mobility of multi-cluster is- 
lands. 



II. COMPUTATIONAL DETAILS 

Diffusion coefficients for clusters can only be obtained 
at the expense of very long MD runs: there exists nu- 
merous possible diffusion paths, and there is therefore 
not a single energy barrier (and prefactor) characteriz- 
ing the dynamics. These systems, further, do not tead, 
themselves readily to accelerated MD algorithms. E3'E2l 
Brute-force simulations - long enough for statistically- 
significant data to be collected - therefore appear to 
be the only avenue. This rules out ab initio methods, 
which can only deal with very small systems (a few tens 
of atoms) over limited timescales (tens of picoseconds at 
best): empirical or semi-empirical potentials must be em- 
ployed. 

As mentioned above, we describe here the interactions 
betweeii—^u particles using the embedded-atom method 
(EAM),l3 an n-body potential with proven ability to 
describe reliably various static and dynamic properties 
of transition and noble metals in either bulk or surface 
configurations .Ej The model is "semi-empirical" in the 
sense that it approaches the total-energy problem from 
a local electron-density viewpoint, but using a functional 
form with parameters fitted to experiment (equilibrium 
lattice constant, sublimation energy, bulk modulus, elas- 



tic constants, etc.). 

The interactions between C atoms are modelled us- 
ing the Tersoff potential,u an empirical n-body poten- 
tial which accounts well for various conformations of car- 
bon. The Tersoff potential for carbon is truncated at 2.10 
A, which turns out to be smaller than the inter-plane 
distance in graphite, 3.35 A. Thus, within this model, 
there are no interactions between neighbouring graphite 
planes. This is of course an approximation, but not a bad 
one since basal planes in graphite are known to interact 
weekly. (This is why it is a good lubricant!). A pleas- 
ant consequence of this is that the substrate can be as- 
sumed to consist of a single and only layer, thus reducing 
formidably the (nevertheless very heavy) computational 
load of the calculations. 

Last — and most problematic — is the Au-C interac- 
tion, for which no simple (empirical or semi-empirical) 
model is to our knowledge available. One way of de- 
termining this would be to fit an ab initio database to 
a proper, manageable functional potential. However, 
since Au-C pairs conform in so many different ways in 
the present problem, this appears to be a hopeless task, 
not worth the effort in view of the other approximations 
we have to live with. We therefore improvised this in- 
teraction a little bit and took it to be of the Lennard- 
Jones form, with a — 2.74 A and e = 0.022 eV, trun- 
cated at 4.50 A. The parameters were determined rather 
loosely fromAiarious two-body models for Ag-C and Pt-C 
interactions .ES Overall, we expect our model to provide 
a qualitatively correct description of the system, realistic 
in that the most important physical characteristics are 
well taken into account. It is however not expected to 
provide a quantitatively precise account of the particu- 
lar system under consideration, but should be relevant 
to several types of metallic clusters which bind weekly to 
graphite. 

We consider here gold nanoclusters comprising 249 
atoms, a size which is. close to that of clusters deposited in 
the experiments.EHfcJ The graphite layer has dimensions 
66.15x63.65 A^ and contains 1500 atoms. Calculations 
were carried out for several temperatures in the range 
400-900 K. It should be noted that a free-standing 24^ 
atom Au cluster melts at about 650 K in this model.El 
This temperature is not affected in a significant man- 
ner by the graphite substrate as the interaction between 
Au and graphite-C atoms are weak. However, the dy- 
namics of the cluster is expected to be different in the 
high-temperature molten state and the low-temperature 
solid state. All simulations were microcanonical, except 
for the initial thermalisation period at each temperature; 
no drift in the temperature was observed. 

Simulations were carried out in most cases using a fully 
dynamical substrate. In two cases, one for the single 
cluster and the other for the dicluster, extremely long 
runs using a static (frozen) substrate were performed: 
it has been found by Deltour et all! that diffusion is 
quantit atively similar on both substrates (however, see 
Section IIIB below). The equations of motion were in- 
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tegrated using the velocity form of the Verlet algorithm 
with a timestep of 1.0 asd 2.5 is for dynamic and static 
substrates, respectively.Ea (Carbon being a light atom, a 
smaller timestep is needed in order to properly describe 
the motion). The dynamic-substrate simulations ran be- 
tween 10 and 14 million timesteps (depending on tem- 
perature), i.e., 10-14 ns. The static-substrate simulation 
for the monocluster, in comparison, ran for a total of 50 
million timesteps, i.e., a very respectable 125 ns = 0.125 
/is; the corresponding dicluster simulation ran for 75 ns. 
All calculations were performed using the program groF, 
a general-purpose MD code for bulk and surfaces devel- 
oped by one of the authors (LJL). 



III. RESULTS 
A. Dynamic-substrate simulations 

We first discuss diffusion on a dynamic substrate, i.e., 
with all parts of the system explicitly dealt with in the 
MD simulations. Fig. |l| gives the (time-averaged) mean- 
square displacements (MSB's) of the cluster's center-of- 
mass at the various temperatures investigated, which 
will be used to calculate the diffusion constant, D — 
limt_+oo r'^{t)/4:t. As indicated above, the simulations ex- 
tend over 10-14 ns, but the MSB's are only shown for a 
maximum correlation time of one ns in order to "ensure" 
statistical reliability. It is evident (e.g., upon comparing 
the results at 700 and 800 K) that the diffusion coeffi- 
cients that can be extracted from these plots will carry 
a sizeable error bar. Nevertheless, it is certainly the case 
that (i) diffusion is very significant and (ii) it increases 
rapidly with temperature. 

There is no evidence from these plots that the MSB's 
obey a non- linear power law behaviour (i.e., that the clus- 
ter undergoes superdiffusion) which could be associated 
with "Levy flights" : the statistical accuracy of the data is 
simply not sufhcient to draw any conclusions. The clus- 
ter does however undergo long jumps during the course 
of its motion. We will return to this point below when 
we discuss diffusion on a frozen substrate. 

In lack of a better description of the long-time be- 
haviour of the diffusion process, we simply assume that 
r'^{t) — * 4:Dt as t gets large. The resulting diffusion coef- 
ficients are plotted in the manner of Arrhenius, i.e., logZ? 
vs l/ksT, in the inset of Fig. |]. If the process were truly 
Arrhenius, all points would fall on a single straight line. 
This is evidently not the case here. Though we could 
probably go ahead and fit the data to a straight line, 
attributing the discrepancies to statistical error, there is 
probably a natural explanation for the "break" that a 
sharp eye can observe between 600 and 700 K: As noted 
above, the free A1U249 cluster melts at about 650 K in 
the EAM modeLQ The presence of the substrate raises 
the melting point, but very little since the interactions 
between the cluster and the graphite surface are small. 



Thus, the cluster is solid at the lowest temperatures (400, 
500 and 600 K), but liquid above (700, 800 and 900 K). 
The statistics are evidently insufficient to allow firm con- 
clusions to be drawn; there nevertheless appears to be 
a discontinuity near the cluster melting point tempera- 
ture, with activation energ ies on either side of about 0.05 
eV. We discuss in Section [II B the implications of these 
findings on the kinetics of growth. 



B. Static-substrate simulations 

The static-substrate simulations, carried out at a sin- 
gle temperature (for the cluster), viz. 500 K, serve many 
purposes: (i) re-assess the equivalence with |-dynamic- 
substrate MB runs reported by Beltour et al.p (ii) pro- 
vide accurate statistics for a proper comparison of the 
diffusive behaviour of mono- and diclusters; (iii) examine 
the possible superdiffusive character of the trajectories. 

We focus, first, on a comparison between static- and 
dynamic-substrate simulations. As can be appreciated 
from the MSB's given in the inset of Fig. ^, there is 
a rather substantial difference between the two calcula- 
tions: for the dynamic substrate at 500 K, the diffusion 
constant is 3.71 x 10~^ cm^/s, while for the frozen sub- 
strate we have 1.09 x 10~^ cm^/s. (This value is actually 
significantly smaller than that for the dynamic substrate 
at 400 K — 100 K lower temperature — viz. 1.70 x 10~^ 
cm^/s). Again, statistical uncertainties cannot be totally 
excluded to account for this discrepancy, but it is difh- 
cult to imagine that it could explain all of the observed 
difference (cf. inset to Fig. |l] for a better appreciation of 
this difference) . The explanation might however be quite 
simple. 

As noted above, the cluster-substrate interactions are 
weak, and this likely plays an important role in determin- 
ing the characteristics of the motion. Visual inspection 
of the X — y paths in the two different situations makes 
it apparent that the motion has a much stronger "stick- 
and-jump" character on the frozen substrate than on the 
dynamic one. On the frozen substrate, further, the tra- 
jectory is more compact on a given timescale. This can 
in fact be characterized in a quantitative roanner by con- 
sidering, following Luedtke and Landmanjj the function 
Prid), which gives the distribution of displacements of 
length d over a timescale of r. The motion is best char- 
acterized using a value of r corresponding to the period 
of vibration of the cluster in a sticking mode (see below). 

The function Pr^d) (normalized to unity) is displayed 
in Fig. H for the dynamic substrate at three different tem- 
peratures (400, 500, and 900 K) and for the static sub- 
strate at 500 K. The value of t was determined from 
the frozen-substrate simulations by simply counting the 
number of oscillations over a given period of time; we 
found T = 20 ps to within about 10%. We note that, for 
the dynamic substrate, the period of oscillations at 400 
K is about 38 ps, while no oscillations can be found at 
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500 and higher temperatures, i.e., the sticking mode is 
absent above 500 K or so. 

The difference between static and dynamic substrates 
is striking: On the dynamic surface, Prid) is a broad 
featureless distribution, which gets broader as tempera- 
ture increases. The maximum of the distribution at low 
temperature lies at about 1.6-1.8 A — roughly the dis- 
tance between equilibrium sites on the graphite surface 
— clearly establishing that the motion proceeds in an 
quasi continuous manner via "sliding hops" to nearest- 
neighbours; the hops get longer as temperature increases. 
On the static substrate, in contrast, a "sticky" vibra- 
tional mode, of amplitude roughly 0.25 A, is clearly vis- 
ible. This is followed by a broad tail which corresponds, 
again, to the sliding jumps that are characteristic of the 
motion on the dynamic substrate. 

Sticking, therefore, is much more likely to take place on 
the static than on the dynamic substrate, thereby con- 
tributing to decrease the average distance traveled by 
the cluster over a given period of time. This conclu- 
sion is however not general: The system under consider- 
ation here is perhaps a bit peculiar in that the cluster- 
substrate interactions are especially weak. (In compari- 
son, Luedtke and Landman's e for the Au-C interaction is 
0.01273 eV, even smaller than our own value.) One may 
conjecture that the vibrations of the surface are enough, 
in such cases, to overcome completely the barrier oppos- 
ing diffusion, which might not be true of systems where 
the interactions are stronger (as in the case, e.g., of Del- 
tour et al.'s simulations, Ref. [7|). 

It also appears that our diffusion data do not cover a 
timescale long enough to warrant firm conclusions to be 
drawn on the possibility that superdifFusion might be tak- 
ing place. This is certainly true, as we have seen above, 
of the dynamic-substrate simulations, which cover "only" 
10-14 nanoseconds, but also of the static-substrate sim- 
ulations (assuming, in view of the above discussion, that 
they are relevant to the problem under study), which ex- 
tend to 125 ns. Certainly, the position of the cluster's 
center-of-mass does exhibit something of a self-similar 
character, as reported by Luedtke and Landmadj, and as 
can be seen in Fig. |4[ Evidently, one cannot trust statis- 
tics here over more than a decade or two in time. One 
might hope that superdiffusion would be more apparent 
in the long-time behaviour of the MSB's. To this effect, 
we plot in Fig. |^ log (Or^ lc>§ ^r a maximum correla- 
tion time of (here) 40 ns;c3 the slope of such a plot is the 
diffusivity exponent 7. The statistical quality of the data 
decreases with correlation time, and becomes clearly in- 
sufficient over 5 ns or so; the large dip at about 15 ns 
can testify. Our best estimate of the slope 7 at "large" 
(more than ~1 ns) correlation times is anywhere between 
0.9 and 1.2, i.e., mild underdiffusion or mild superdiffu- 
sion... or no superdiffusion at all! This is consistent 
with the value reported by Luedtke and Landman, who 
find 7 = 1.1 based on an analysis of sticking and sliding 
times. One point worth mentioning is that the velocity- 
autocorrelation function for adatom diffusion in the in- 



termediate and high-friction regimes has been shown to 
follow a power-law behaviour at intermediate timesj-the 
exponential dependence resumes at very long timesES 

C. Diclusters 

The morphology of films grown by cluster deposition 
depends critically on the coefficient of diffusion of mono- 
clusters, as we have just seen, but also, because clusters 
aggregate, on the coefficient of diffusion of multiclusters. 
From simple geometric arguments, it might be argued 
that the rate of diffusion should scale as iV~^/'^, where 
N is the number of atomsjn the cluster, as was in fact 
observed by Deltour et alB for Lennard- Jones clusters. 
However, it can be expected that the morphology of the 
films depends, as well, on the shape of the multi-clusters 
following the aggregation of monoclusters, i.e., on the 
kinetics of coalescence. .-j 

In a previous publicationjj we examined the coales- 
cence of gold nanoclusters in vacuum and found it to 
be much slower than predicted by macroscopic theories. 
This state of affairs can be attributed to the presence of 
facets and edges which constitute barriers to the trapsr 
port of particles required for coalescence to take placeEa. 
The "neck" between two particles was however found to 
form very rapidly. We conjectured that these conclu- 
sions would apply equally well to the particular case of 
gold nanoclusters on graphite since the gold-graphite in- 
teractions are weak. 

We have verified this in the context of the present work: 
indeed, coalescence is little affected by the presence of the 
substrate, as demonstrated in Fig. |[ We considered both 
a free-standing and a supported pair of 249-atom gold 
clusters. Starting at very low temperature (50 K), tem- 
perature was slowly and progressively (stepwise) raised 
to 600 K. (As noted above, the 249-atom gold cluster 
melts at about 650 K in this EAM model and we there- 
fore did not go beyond this point). We plot, in Fig. ||, the 
evolution with time-temperature of the three moments of 
inertia of the dicluster. Since the cluster can rotate, the 
moments of inertia provide a more useful measure ofJhe 
shape of the object than, e.g., the radii of gyrationjj A 
side view of the dicluster at 200 K, i.e., after the neck 
between the two monoclusters has formed completely, is 
shown in Fig. ||. It is evident that the dicluster does 
not wet the surface, and therefore the substrate plays a 
relatively minor role in the coalescence process. 

As can be seen in Fig. ^, the behaviour of the free- 
standing and supported diclusters are almost identical, 
except for the initial phase of coalescence: the sup- 
ported cluster forms a neck much more rapidly than the 
free-standing cluster, presumably because the substrate 
offers, through some thermostatic effect, an additional 
route via which coalescence (by plastic deformation) can 
be mediated; it is conceivable also that the substrate 
"forces" the atomic planes from the two clusters to align. 
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We have not explored these questions further; it remains 
that the end points of the two coalescence runs are iden- 
tical within statistical uncertainty. Thus, again, coales- 
cence is hampered by the presence of facets and edges; the 
timescale for complete coalescence is much longer than 
predicted by continuous theories. The shape of islands 
on the graphite surface will be strongly affected, and it 
is also expected that the rate of diffusion will be affected 
(since it is determined by the contact area between sub- 
strate and cluster). 

The MSD of the dicluster (after proper equilibration at 
500 K) is displayed in the inset of Fig. ||. As mentioned 
earlier, this was calculated from a static-substrate run 
covering 75 ns. The same limitations as noted above for 
the monocluster should therefore hold in the present case. 
It is a very remarkable (and perhaps even surprising) 
result that the rate of diffusion of the dicluster is quite 
comparable to that of the monocluster, inasmuch as the 
frozen-substrate simulations are concerned. (We expect 
the diffusion constants on the dynamic substrate to be 
different — and larger — but in a proportion that would 
be quite comparable to that found here). The value of 
D = 1.38 X 10~^ cm^/s we obtain for the dicluster is in 
fact a bit larger than that for the monocluster (1.09 x 
10^^ cm^/s). The difference is probably not meaningful; 
what is meaningful, however, is that the the mono- and 
the dicluster have comparable coefficients of diffusion; 
this has p rofoun d implications on growth, as we discuss 
in Section III D| , below. 

The function P{d) for the dicluster at 500 K is dis- 
played in Fig. |3[ here we estimate that r = 40 ps (vs 
about 20 ps for the monocluster). The distribution is 
quite similar to that found for the single cluster on the 
frozen substrate, though broader and shifted to slighlty 
larger displacements. This last result is likely due to 
the fact that, being larger, the dicluster is not as easily 
able to accomodate itself with the substrate as the mo- 
nocluster; in this sense, it is more loosely bound to the 
substrate. 



D. Comparison with experimental results 

Deposition of gold clusters on graph ite experiments 
were carried out in Lyon recently.ll3llj Several models 
have been proposed to extract the microscopic cluster di£i 
fusion coefficients from the measured island densities.tHI 
Of course, in order to provide a meaningful interpretation 
of the data, the models must take into account the precise 
conditions in which the experiments are performed. In 
Lyon, for instance, the flux of clusters is chopped, rather 
than continuous, and this-affiacts the kinetics of diffusion 
and growth considerablyJ23'Ell Previous estimates of the 
rates of diffusion of Au on graphite, which overlooked 
this important detail, are therefore in error. In Ref. |l^, 
a diffusion coefficient of lO"'^ cm^/s at 400 K is given; 
for a discussion, see Ref. 12l The "correct" number, in- 



cluding flux chopping, would be 1.0 cm^/s if monoclus- 
ters only were assumed to he mobile. However, as we 
have seen above, cluster dimers diffuse at a rate which 
is quite comparable to that for monoclusters, suggesting 
that larger clusters would diffuse |as well. The Lennard- 
Jones simulations of Dcltour et aln indicate that the rate 
of diffusion of compact A^-atom clusters scales roughly as 
the inverse of the contact area between the cluster and 
the substrate: Dn — DiN^^^"^. (Compact clusters are 
expected to form through aggregation and coalescence; 
see Ref. |l2|). 

Experimentally, however, it is almost impossible to de- 
termine whether or not multiclusters do diffuse, and at 
which rate. In view of this, and the expected importance 
of multicluster mobility on growth, we have carried out 
a series of kinetic Monte Carlo (KMC) simulations in or- 
der to estimate the largest island which must be allowed 
to diffuse in order to account for the experimentally- 
observed gold island density on graphite Mr^PO K, viz. 
4 X 10^ islands/cm^, or 1.1 x 10^^ per site.yy To do so, 
we assume that the diffusion constant for monoclusters 
found in the present simulations is correct, and that the 
rate of diffusion of A^-clusters scales according to the law 
given above. All other parameters (incident cluster flux, 
temperature, chopping rate, etc.) are fixed by experi- 
ment. 

Figure ^ shows the results of the KMC simulations: 
we plot here the island density that would be observed if 
the largest mobile island were of size A^max- The com- 
putational load increases very rapidly with A^max E^nd 
we therefore only considered islands of sizes less than 
or equal to 35. The data points follow very closely a 
power-law relation and we can thus extrapolate to larger 
values of A^max, i.e., smaller island densities. We find in 
this way that islands up to a maximum size of about 100 
mono-clusters must be mobile in order to account for the 
observed island density of 1.1 x lO"'^ per site. In what 
follows, we discuss in more detail the connection of this 
observation with experiment. 

We first note 

that, in the gold-on-graphite experiments,c3 large islands 
form which are "partially ramified" , in the sense that the 
branch width is much larger than the size of the deposited 
clusters, each branch being formed by the coalescence 
of up to 200 monoclusters. In contrast, for antirnaiiy. 
cluster deposition on graphite at room temperatureji3'E3 
the islands are fully ramified, i.e., have a branch width 
identical to the diameter of the monoclusters; this estab- 
lishes unambiguously that cluster coalescence is not tak- 
ing place in this case. It has been shown, further, that 
the mobility of the islands is negligible in antimony.1121 
Our results suggest, therefore, when taken together with 
the work of Deltour et a/.JZI that compact islands, which 
form through diffusion and coalescence, are mobile ac- 
cording to a N~'^/^ law. In contrast, ramified islands, 
which form when coalescence does not take place, have 
much reduced mobility - certainly much less than would 
be expected from a N^^^^ law. A^max, therefore, signals 
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the crossover point between the two mobility regimes or, 
equivalently, the multicluster size at which the morphol- 
ogy of the islands crosses over from compact to ramified 
(or vice- versa). The physical reasons underlying the rela- 
tion between mobility and morphology are not clear, but 
there appears to be no other ways to interpret the ex- 
perimental results. This problem clearly deserves further 
studies. 

To summarize this section, the mobility of large is- 
lands is evidently a necessary ingredient to account for 
the experimentally observed island density. Our simula- 
tions suggest that these islands can be as large as 100 
monoclusters; while this is consistent with experiment, 
the exact value, as well as the precise dependence of the 
diffusion rate on size, cannot at present be estimated. 

IV. CONCLUDING REMARKS 

Cluster-deposition techniques are of great potential 
interest for assembling materials with specific, tailor- 
made applications. Yet, the fabrication process depends 
critically on the possibility for the clusters to diffuse 
on the surface in order to settle in appropriate posi- 
tions, thus forming self-organized structures, or to aggre- 
gate/coalesce with other clusters in order to form larger- 
scale structures and eventually continuous layers. In this 
article, we have demonstrated, using molecular-dynamics 
simulations with realistic interatomic potentials, that the 
diffusion of large metallic clusters on graphite can take 
place at a pace which is quite comparable to that for 
single adatoms. We have also established that the rate 
of diffusion of cluster dimers can be very sizeable, com- 
parable in fact to that for monoclusters. An extremely 
important consequence of this is that islands formed by 
the aggregation of clusters are also expected to be mo- 
bile. Using kinetic Monte Carlo simulations and assum- 
ing a proper scaling law for the dependence on size of 
the diffusivity of large clusters, we estimate that islands 
containing as much as 2-5 000 atoms (100 monoclusters) 
are expected to undergo diffusion at a significant rate 
on graphite surfaces. These findings have profound con- 
sequences for the morphology of cluster-assembled thin 
films. 
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FIG. 2. Main figure: Log-log plot of the time-averaged 
mean-square displacements for the cluster's center-of-mass 
on the static substrate at 500 K. The three curves corre- 
spond to different estimates: using the full extent of the 
!=«iie(frfli"«&rve); only the first half (dashes); only the sec- 
ond half (dots). The difference between these curves gives a 
measure of the error on the estimated diffusion coefficient. In- 
set: Time-averaged mean-square displacements for the mon- 
ocluster on a static substrate (full line), the monocluster on 
a dynamic substrate (dashes) and the dicluster on a static 
substrate (dots). 
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FIG. 1. Main figure: Time-averaged mean-square dis- 
placements for the cluster's center-of-mass at the various tem- 
peratures investigated, namely, from bottom to top, 400, 500, 
600, 700, 800 and 900 K. (The 700 and 800 K curves are in- 
Veit'Sd. fopefl > 0.7 ns). Inset: Arrhenius plot of the diffusion 
coefficient. The open square at 500 K {l/ksT = 23.2 eV~^) 
is the result for the frozen substrate. 
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FIG. 3. The function Prid), which gives the distribution 
of displacements of length d over a timescale of r, at three 
different temperatures, as indicated. At 500 K, the three 
curves correspond to the monocluster on a dynamic substrate 
(ftilPHHSJpf'^ftie monocluster on a static substrate (dashes) and 
the dicluster on a static substrate (dots). 
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FIG. 4. The x position of the cluster's center of mass 
at 500 K on a frozen substrate for three different timescales, 
showing the apparently self-similar character of the trajectory. 



FIG. 5. Evolution with time-temperature of the three mo- 
ments of inertia (full lines) of the dicluster, both free-standing 
and supported on graphite; the correspondence between tem- 
perature and time is indicated by the dashed, stepwise curve. 




FIG. 6. Ball-and stick model of the gold dicluster on the 
graphite substrate at 200 K, after the neck between the two 
monoclusters has formed completely. The two monoclusters 
are colored differently for ease of visualisation. 



n ' \ ' \ ' r 

Free standing 



0.4 



Lewis et al, Figure 5 




2 3 4 5 6 
t (ns) 



u 



m 
C 

T3 

CO 



tlfl 
o 




log N„ 



FIG. 7. Predicted island density as a function of the size 
A'^max of the largest multi-cluster island which is allowed to 
diffuse. The calculations were carried out using kinetic Monte 
KSfcf'si'SiuIations as discussed in the text. The full line is a 
linear fit to the data points and the dashed line indicates the 
experimental density. 
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